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An analysis of adiabatic perturbations of a perfect fluid is performed to first-order about a general 
FLRW background using the 1 + 3 covariant and gauge-invariant formalism. The analog of the 
Mukhanov-Sasaki variable and the canonical variables needed to quantise respectively the scalar 
and tensor perturbations in a general FLRW background space-time are identified. The dynamics 
of the vector perturbations is also discussed. 



I. INTRODUCTION 

Cosmological perturbation theory about a given back- 
ground metric was initiated by Lifshitz [1] and extended 
notably by Bardeen [2J, Mukhanov, Feldmann & Bran- 
denberger [3J, and Kodama & Sasaki jj]. This is the 
standard approach to performing a perturbation analy- 
sis of general relativity in order to describe the 'realistic' 
dynamics of cosmological models. It has the disadvan- 
tage that truly physical results can be obtained only af- 
ter completely specifying the correspondence between the 
'real' perturbed space-time and a 'fictitious' unperturbed 
background space-time, which is chosen to be highly sym- 
metric. Such a correspondence is not uniquely defined 
and changes under a gauge transformation. The degrees 
of freedom in the definition of the correspondence, or the 
gauge freedom, leave unphysical gauge modes in the dy- 
namical equations describing the evolution of the pertur- 
bations [5 a . This gauge problem [6] |7] is inherent to such 
a 'background-based' perturbation approach. Indeed, the 
metric, and consequently the Einstein equations, have 10 
degrees of freedom whereas the dynamics is determined 
by 6 parameters only. Therefore, the 4 remaining degrees 
of freedom are directly related to the gauge. 

To solve this problem, Bardeen [2] determined a set of 
gauge-invariant quantities to describe the perturbations 
and derived their dynamical equations. These quantities 
are mathematically well defined but do not have a trans- 
parent geometrical meaning since they are defined with 
respect to a particular coordinate system [5], and their 
physical meaning is obscure 

An alternative way to circumvent the gauge problem 
is to follow the 1 + 3 covariant approach, which was de- 
veloped by Hawking [TU] and extended by Olsen [TT] and 
Ellis & Bruni [12]. The aim of this approach is to study 
the dynamics of real cosmological fluid models in a phys- 
ically transparent manner. This formalism relies on co- 
variantly defined variables, which are gauge-invariant by 
construction [5], thus simplifying the methodology and 
clarifying the physical interpretation of the models. It 
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also allows the metric to be arbitrary. This approach 
admits a covariant and gauge-invariant linearisation that 
allows a perturbation analysis to be performed in a direct 
manner |13j . 

The basic philosophy of the perturbation theory based 
on the 1 + 3 covariant formalism is different from the 
'background-based' perturbation theory [13] . Instead of 
starting with a background space-time and then perturb- 
ing it, the approach begins with an inhomogeneous and 
anisotropic (real) space-time, which reduces to the back- 
ground space-time on large scales. In our case, we will 
take the background to be a homogeneous and isotropic 
FLRW space-time, although this formalism allows in 
principle for more complicated backgrounds. Therefore, 
the 'real' space-time has been appropriately called an al- 
most FLRW space-time [§]. In this perturbation the- 
ory, the approximation takes place by neglecting higher- 
order terms in the exact equations when the values of the 
kinematic and dynamic variables are close to those they 
would take in the background FLRW space-time. The 
analysis is performed in the real space-time and the dy- 
namical equations are subsequently linearised. The back- 
ground solution is simply the zero-order approximation 
of the exact solution of the dynamical equations. 

Although in the 1 + 3 covariant formalism, the scalar, 
vector and tensor perturbations are handled in a unified 
way [13], they decouple to first-order and can be stud- 
ied independently [3J. For a perfect fluid in a spatially- 
curved case, the Bardeen variable for the scalar perturba- 
tions in the 1 + 3 covariant approach was first identified 
by Woszczyna & Kulak [TS] and the curvature pertur- 
bation by Bruni, Dunsby & Ellis 5 . The first attempt 
to establish an explicit relation between the perturbation 
variables in the 1 + 3 covariant formalism and the cor- 
responding variables in the 'background-based' approach 
was made by Goode [Trj] , 

In the 1 + 3 formalism, the time derivative of a phys- 
ical quantity is defined usually as the projection of the 
covariant derivative of the quantity on the worldline as 
outlined in Sec. [Tl] but this is not the only way to define 
covariantly a time derivative as Thiffeault showed in |17j . 
As Langlois & Vernizzi suggested [TB], the Lie derivative 
along the worldline of a fluid element is another possi- 
ble definition. However, the Lie derivative of a scalar 
field along the worldline is identical to time derivative of 
the scalar field along the worldline. Since, in this pub- 



2 



lication, we aim to identify scalar quantities, which are 
respectively scalar perturbations and the scalar ampli- 
tude of vector and tensor perturbations, we will define 
the time derivative of physical quantities to be the pro- 
jection of the covariant derivative of the quantity on the 
worldlinc. 

It is worth mentioning that, recently, Pitrou & 
Uzan [JJ2], working in a spatially-flat case, used Lie 
derivatives to recast the dynamical equations for the 
scalar and tensor perturbations in the 1 + 3 covariant for- 
malism in order to identify perturbation variables which 
are similar to the Sasaki-Mukhanov and canonical vari- 
ables. In particular, they sought to identify the scalar 
and tensor variables that map to the Mukhanov-Sasaki 
and to the tensor canonical variables when considering 
a spatially-flat almost FLRW universe. In the current 
work, we seek the extension to general FLRW universes. 
However, there are some further differences with the work 
of Pitrou & Uzan. Firstly, in the scalar case, the pertur- 
bation variable v a that Pitrou & Uzan obtained is not a 
scalar but a covector. We will seek a scalar here. Sec- 
ondly, for tensor perturbations in the spatially-flat case, 
Pitrou & Uzan obtained a wave equation in terms of 
Lie derivatives of tensor fields. In the current work, we 
will extend their analysis by seeking to identify canoni- 
cal variables, which correspond to the scalar amplitudes 
of the tensor perturbations. These scalar amplitudes are 
the canonical variables defined by Grishchuk 20J. We 
will be using the name 'canonical' to refer to these vari- 
ables henceforth in this publication. Finally, Pitrou & 
Uzan did not consider a perfect fluid analysis as we carry 
out here, but instead restricted their analysis to a single 
scalar field. 

Thus the aim of this publication is to identify, in the 
1 + 3 covariant approach, for adiabatic perturbations of 
a perfect fluid in spatially-curved FLRW models, the 
analog of the Mukhanov-Sasaki variable [3J 0] and the 
canonical variables |20j needed to quantise respectively 
the scalar and tensor perturbations. The quantisation 
of fields in a spatially-curved space-time lies outside the 
scope of this paper. However, this topic has been ad- 
dressed in detail by Birrell & Davis [21] and Fulling [22]. 
Our approach broadens the scalar perturbations analysis 
performed by Woszczyna & Kulak |15j , Bruni, Dunsby & 
Ellis [5] and Lyth & Woszczyna [53] in a spatially-curved 
case. 

The structure of this publication is as follows. In 
Sec. [Hj we give a concise description of the kinematics 
in the 1 + 3 covariant approach. In Sec. |III[ we establish 
the first-order dynamical equations for adiabatic pertur- 
bations of a perfect fluid in the 1+3 covariant approach. 



In Sec. IV we identify the analog of the Mukhanov-Sasaki 
variable for the scalar perturbations. In Sec.[V] we deter- 
mine the dynamics of the scalar amplitude of the vector 



perturbations. In Sec. VI we finally identify the canoni- 



cal variables for the tensor perturbations. 

For convenience, we follow Hawking [TO] and Ellis [T5] 
by adopting the (—,+,+,+) -signature for the metric. 



This choice of signature is particularly appropriate for 
the 1 + 3 covariant formalism. Indeed, the dynamical 
equations are projected on the local spatial hypersur- 
faces, which are positively defined for our signature con- 
vention. The correspondence for dynamical quantities 
expressed in terms of the opposite signature can be found 



II. KINEMATICS IN THE 1+3 COVARIANT 
FORMALISM 



We will briefly outline the basics of the 1 + 3 covariant 
formalism describing the fluid kinematics. To introduce 
this formalism, we follow Ellis & van Elst's approach |25j . 
The approach is based on the 1 + 3 decomposition of 
geometric quantities with respect to a fundamental 4- 
vclocity u a that uniquely determines the frame and the 
worldline of every infinitesimal volume element of fluid, 



dx a 



u a u 



-1 



(1) 



where x a are arbitrary cosmic coordinates, and r is the 
proper time measured along the worldlines. In the con- 
text of a general cosmological model, we require that the 
4-velocity be chosen in a physical manner such that in 
the FLRW limit the dipolc of the cosmic microwave back- 
ground radiation vanishes. This condition is necessary to 
ensure the gauge-invariance of the approach. 

The 4-velocity u a defines locally and in a unique fash- 
ion two projection tensors, 



Uab = -u a u b , 

hab = 9ab + U a U b , 

which satisfy the following properties, 

U a c U c b = U% , U\ = 1 , U ab u b = u a 
h\h c b = h\ , h\ = 3 , h ab u b = . 



(2) 



(3) 



The first projects parallel to the 4-velocity vector u a , 
and the second determines the metric properties of the 
instantaneous rest frame of the fluid. It is useful to intro- 
duce the totally antisymmetric Levi-Civita tensor £ abc , 
which is defined in the rest frame of the fluid and thus 
satisfies, 



£ab c u c = 



(4) 



Moreover, we define two projected covariant deriva- 
tives which are the time projected covariant derivative 
along the worldline (denoted ' ) and the spatially pro- 
jected covariant derivative (denoted D a ). For any quan- 
tity Q a - - b , these are respectively defined as 



Q a - b 

D C Q° 



u c V c Q a - 
= h f r h a r . 



b ■ ■ 



,/... 



(5) 
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Furthermore, the kinematics and the dynamics are de- 
termined by projected tensors that are orthogonal to u a 
on every index. The angle brackets are used to denote 
respectively orthogonal projections of vectors V a and the 
orthogonally projected symmetric trace-free part (PSTF) 
of rank-2 tensors T ab according to, 



V {a) = h a b V b , 



T 



cd 



(6) 



For convenience, the angle brackets are also used to de- 
note the orthogonal projections of covariant time deriva- 
tives of vectors and tensors along the worldline u a as 
follows, 



V {a) = h a b V b 



f{ab) = ^ h {a ch b) d _l h ab hci]T 



(7) 



Note that, in general, the time derivative of vectors and 
tensors does not commute with the projection of these 
quantities on the spatial hypersurfaces according to, 



(8) 



The projection of the covariant time derivative of a quan- 
tity Q a " b on the spatial hypersurfaces is defined as, 

W(Q a - b ..J =h a c ...h\...u e V e Q c - d „, . (9) 

It is also useful to define the projected covariant curl as, 

cur\Q a ... b = s cd{a D c Q d ^ b) . (10) 

Information relating to the kinematics is contained in 
the covariant derivative of u a which can be split into 
irreducible parts, defined by their symmetry properties, 



V a ti & = ~u a a b + D a u b 



-u a a b + \<dh ab + a ab + LO ab 



(11) 



where 



• a a ee u b V b u a is the relativistic acceleration vector, 
representing the degree to which matter moves un- 
der forces other than gravity. 

• ee D a u a is the scalar describing the volume rate 
of expansion of the fluid (with H = §6 the Hubble 
parameter) . 

• a ab ee D( a u b } is the trace-free rate-of-shear tensor 
describing the rate of distortion of the fluid flow. 

• io ab ee D[ a ii b ] is the antisymmetric vorticity tensor 
describing the rotation of the fluid relative to a non- 
rotating frame. 



These kinematical quantities have the following proper- 
ties, 



a a u a = , 

CTafcU 6 = , a ab — <7( tt j,) 

uj ab u b = , uj ab = uj [ab] 



(j 



(12) 



Note that in presence of vorticity the spatially pro- 
jected covariant derivatives do not commute, which im- 
plies that any scalar field S has to satisfy the non- 
commutation relation, 



D[ a D b ]S — uj ab S . 



(13) 



It is useful to introduce a vorticity pseudovector iv a , 
which is the dual of the vorticity tensor uj bc and is defined 

as, 



, a — 1 -ahc 

lu = uj bc , 



and a vorticity scalar given by, 

^ = M) 1/2 



(14) 



(15) 



To build a covariant linear perturbation theory, we will 
linearise the quantities of the 1 + 3 covariant formalism 
about an FLRW background space-time which is, by def- 
inition, homogeneous and isotropic. Gauge invariance of 
the perturbed quantities is guaranteed by the Gauge In- 
variance Lemma, which states that if a quantity vanishes 
in the background space-time, then it is gauge invariant 
at first-order The homogeneity of the background 

space-time implies that the acceleration a a is a first-order 
variable since it vanishes in the background space-time 
(i.e. at zeroth order). Similarly, the isotropy of the back- 
ground space-time implies that rate of shear a ab is also 
a first-order variable. Finally, the existence of hyper- 
surfaces orthogonal to the worldline in the background 
space-time implies that vorticity tensor ui ab and the vor- 
ticity covector uj a are first-order variables. Therefore, the 
only zero-order kinematic quantity is the rate of expan- 
sion O. 

In the background space-time, which means to zero- 
order in the dynamical variables, the covariant derivative 
of the worldline thus becomes, 



y a u b 



^Qh ab 



(16) 



It is also useful to define (up to some constant factor) 
a zero-order scale factor R such that, 



e = 3H 



R 
R 



(17) 



where H is the cosmic Hubble scale factor. 

Finally, it is also convenient to introduce a conformal 
time variable f satisfying the differential relation, 



dr 



dr 
~R 



(18) 
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which implies that the derivatives with respect to cosmic 
time r and conformal time f of a quantity Q a ' b ... are 
related by, 



RQ a - b 



(19) 



where a prime denotes a derivative with respect to con- 
formal time f . The conformal Hubble scale factor H is 
related to the cosmic scale factor H by H = RH. 



III. FIRST-ORDER ADIABATIC DYNAMICS 
OF A PERFECT FLUID ABOUT AN FLRW 
BACKGROUND 

We will now use the 1 + 3 covariant formalism, out- 
lined in Sec. |TT[ to describe the adiabatic dynamics of 
a perfect fluid to first-order in the dynamical variables, 
which means neglecting second- and higher-order prod- 
ucts of first-order dynamical variables. We then per- 
form a perturbation analysis of such a fluid about an 
FLRW background, which will be used in Sec. |IV||VI| 
to describe scalar, vector and tensor perturbations re- 
spectively. Thus, we require the cosmological fluid to be 
highly symmetric on large scales but allow for generic 
inhomogeneities on small scales. 

The dynamics of a perfect fluid is described by the 
Einstein field equations, which read, 



Rab ndabR- — K T a , 



(20) 



where R ab and 1Z are respectively the Ricci tensor and 
scalar. The dynamical model is fully determined by the 
matter content and the curvature. The matter content 
is described by the stress-energy momentum tensor T ab . 
For a perfect fluid, using the 1 + 3 formalism, it can be 
recast as, 



Tab = PU a U b + ph ab 



where p is the energy density and p the pressure of 
the fluid. We assume the fluid to be a specific linear 
barotropic fluid so that it satisfies the equation-of-state, 



p = wp 



(22) 



where w is the equation-of-state parameter. The energy 
density p and the pressure p are zero-order variables that 
do not vanish on the background. 

For an adiabatic flow, the speed of sound is defined as 



dp 
dp 



(23) 



and the time derivative of the equation-of-state parame- 
ter satisfies, 



w 



-Q{c 2 s -w){l + w) , 



(24) 



and is derived using the energy conservation equa- 



tion (59 1. Note that from (24), it follows that c 



if w = 0. 



IV 



All the information related to the curvature is encoded 
in the Riemann tensor which can be decomposed as |10j . 



R ab cd = C ab c d + 2S^ a [ c R b ^ d ] — ^lZ8 a [ c 5 b ( q 



(25) 



where C ab c d is the Weyl tensor constructed to be the 
trace-free part of the Riemann tensor. 

By analogy to classical electrodynamics, the Weyl ten- 
sor itself can be split, relative to the worldline u a , into 
an 'electric' and a 'magnetic' part |10j according to, 



E a b — C aco dU C U d , 

Hab = C ac bdU C U d — h£adeC de 



bcU 



(26) 
(27) 



where *C a bcd is the dual of the Weyl tensor. Their prop- 
erties follow directly from the symmetries of the Weyl 
tensor, 



E a bU b — , E ab — E( a ty , E a a — , 

HabU b = , H ab = H( a b) , H a a = . 



(28) 



These electric and magnetic parts of the Weyl tensor rep- 
resent the 'free gravitational field', enabling gravitational 
action at a distance and describing tidal forces and grav- 
itational waves. The Weyl tensor vanishes on a confor- 
mally flat background space-time such as FLRW models. 
Thus, the electric E ab and magnetic part H ab of the Weyl 
tensor are first-order variables. 

The Ricci tensor R ab is simply obtained by substitut- 
ing the expression (21 1 for the stress energy momentum 



tensor into the Einstein field equations (|20j), 

Rab = § (p + 3p) u a u b + f (p - p) h ab . (29) 
The Riemann tensor R a bcd can now be recast in terms 



of the Ricci tensor ( 29 1 , the electric ( 26 1 and mag- 



(21) 

netic (p7| parts of the Weyl tensor according to the de 



composition ( 25 1 in the following way, 



R ab cd ^ K (p+3p)u [a u lc h b ^d] + lnph a [c h b d] 

+ ^ a u [c E^ d ]+^ a [c E b \ l] (30) 
+ 2e abe u [c H d]e + 2e cde ^ a H b ^ . 

The Riemann tensor restricted to the spatial hypersur- 
face, ^R ab c d, is related to the Riemann tensor defined 
on the whole space-time, R ab c d, by 

^R ab cd = h a e h b f hS c h h d R ef gh - 2v a lc v b d] . (31) 

where the tensor v ab is defined as, 

v ab ee D a u b . (32) 



For a perfect fluid, using the decomposition (30 1, the 
spatial Riemann tensor (31) becomes, 



^R ab cd - lnph a [c h b d ] + 4h a [c E b d] - 2v a [c v b d] . (33) 
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To zero-order, the decomposition of the spatially pro- 



jected Riemann tensor (31 I reduces to, 



yield the spatial Ricci tensor ^R a b and scalar re- 
spectively, which to first-order are given by, 



(3) nab 2-ff a b 

H cd ~ft2 h \ ch d \ 



(34) 



{3) R i, 

±l ab 



ic = ^e 2 



2np 



\&{<j ab -uj ab ) + E ab , (44) 

(45) 



where if is the curvature parameter, which is related to 
the Gaussian curvature K by, 



To zero-order, using the expression for the Gaussian cur- 
vature (|35| , the contractions of the spatial Ricci identi- 



6K 



(35) 



ties (|44[) and (45) reduce respectively to, 



where is the spatial curvature scalar, which is ob- 
tained by twice contracting the spatially projected Rie- 
mann tensor (311. 



^Rab — § ^TZhab 



K 
R? 



(46) 
(47) 



In general, there are three sets of dynamical equations 
for a perfect fluid. These sets are derived, respectively, 
from the Ricci identities, the Bianchi identities, once- and 
twice-contracted. We present now each set in turn and 
expand the dynamical equations to first-order about an 
FLRW background space-time. 



A. Ricci identities 

The first set of dynamical equations arises from the 
Ricci identities. These identities can firstly be applied to 
the whole space-time and secondly to the spatial hyper- 
surfaces according to, 



Note that the above expression (47) is the Friedmann 



equation. It is useful to recast the Friedmann (47 1 and 



Raychaudhuri (38 1 equations in terms of conformal time. 



To zero-order, these equations are respectively given by, 



V[ a V b ]U c = \R a bcdU d , 



D, n D H v r . = h (3) R ab cdV d 



(36) 
(37) 



n 2 = ipR 2 - k , 

H' = -§pi? 2 (l + 3w) . 

Finally, it is convenient to recast them as, 

H 2 -H' + K= 1pR 2 (l + w) , 
W = -f (l + 3w) (H 2 + K) . 

B. Once-contracted Bianchi identities 



The second and third set of dynamical equations are 
contained in the Bianchi identities. The Riemann tensor 
satisfies the Bianchi identities as follows, 



(48) 
(49) 

(50) 
(51) 



where the spatial vectors v a are orthogonal to the world- 
line, i.e. v a u n = 0. 



V [e R ab] rd = 



(52) 



The information contained in the Ricci identities ( 36 ) 



tions ( 20 1 into t 



re Bianchi identities ( 52 1 and contracting 



By substituting the expression for the Riemann tensor 
(37} can be extracted by projecting them on different hy- decomposition and the effective Einstein field equa 
persurfaces using the decomposition of the corresponding 
Riemann tensors ( 30 1-( 33 ) and following the same proce- 
dure as in 24J. 

To first-order, the Ricci identities applied to the whole 
space-time (36) yield three propagation equations, which 



two indices (d = e), the once-contracted Bianchi identi- 
ties are found to be, 



V d C ab 



are respectively the Raychaudhuri equation, the rate of 
shear propagation equation and the vorticity propagation 
equation, 



cd 



V [a R b] c +lS c [a W b] TZ^0 . (53) 



±e 2 



f (p + Bp) + D b a b , 
\ curl a a , 

& (ab) = — 1@ &ab — E ab + D ( a&) , 
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and three constraint equations, 



D a w n 



, 



D b a n 



|-D a 6 - curlw a , 
H nh = curlcr„(, - Di 



'{a^b) 



(38) 
(39) 
(40) 



(41) 
(42) 
(43) 



In a similar manner to the Ricci identities, the informa- 
tion stored in the once-contracted Bianchi identities has 
to be projected along the worldlines u a and on the spatial 
hypersurfaces h a b . 

To first order, the once-contracted Bianchi identi- 



ties (53 1 yield two propagation equations, which are re- 



spectively the electric and magnetic propagation equa- 
tions, 



E 



(ab) 



= -QE ab + curl H ab ~ f (P + P) (Tab , (54) 



H, 



(ab) 



= -OH ab - curl E ab 



The Ricci identities applied to the spatial hypersur- 
faces express the spatial curvature. Their contractions 



and two constraint equations, 

D b E ab = 1D a p, 
D b H ab = -~n(p + p) ui a 



(55) 

(56) 
(57) 
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C. Twice-contracted Bianchi identities 

The third set of equations is given by the twice- 
contracted Bianchi identities which represent the conser- 
vation of the effective stress energy momentum tensor. 
They are obtained by performing a second contraction 
(6 = c) on the once-contracted Bianchi identities ( 53 ) , 



v & (R ab + y ab n) 



KV b T ab 







To first order, the twice-contracted Bianchi identi- 



ties (58) yield one propagation equations, which is the 



energy conservation equation, 



p = -e (p+ P ) , 



and one constraint equation, which is the momentum 
conservation equation, 



D a p = ~a a (p + p) 



(60) 



IV. SCALAR PERTURBATIONS 

A. Bardeen equation 

Physically, scalar perturbations represent spatial varia- 
tions of zero-order scalar quantities. Thus, spatial Lapla- 
cians of zero-order scalars are natural candidates to de- 
scribe such perturbations. Since we are interested in de- 
termining the time evolution of scalar perturbations in a 
comoving frame, we choose comoving spatial Laplacians 
of zero-order scalars as scalar perturbations variables. 
For the linearised dynamics of a perfect fluid about a 
homogeneous and isotropic background, there are only 
four zero-order scalars, the energy density p, the pres- 
sure p, the expansion rate and the Gaussian curvature 
K. For an adiabatic flow, the pressure is a function of the 
energy density (22 1. Thus, to describe the dynamics of 



adiabatic scalar perturbations, we define three scalar per- 
turbation variables, which are respectively the comoving 
spatial Laplacian of the energy density the comoving 
spatial Laplacian of the expansion rate ^ and the comov- 
ing spatial Laplacian of the Gaussian curvature x, 



$ = RnAp , 
* = RAQ , 



(61) 
(62) 
(63) 



where the comoving spatial Laplacian A is related to the 
cosmic spatial Laplacian D 2 by, 



A = R 2 D 2 = R D a D„ 



(64) 



and the factor of a half in ( 63 1 is included to be consistent 



with the usual definition of the corresponding variable 
in the 'background-based' approach [3]. Similar defini- 
tions for the comoving spatial Laplacian of the energy 



density and the comoving spatial Laplacian of the Gaus- 
sian curvature were used respectively by Woszczyna & 
Kulak 15j and Bruni, Dunsby & Ellis [5]. For a van- 
ishing background curvature, it is worth mentioning that 
the curvature perturbation does not vanish, since it is a 
first-order variable. The dynamics of the scalar pertur- 
bations is obtained by taking the comoving spatial Lapla- 
cian of the scalar propagation equations ( 59 1 , ( 38 1 , ( 45 1 



(5g) and the spatial gradient of the constraint (60 1 in order 



to express the constraint in terms of a comoving spa- 
tial Laplacian of a zero-order scalar. The dynamical 
equations of the scalar perturbations are respectively the 
comoving spatial Laplacian of the energy conservation 
equation ( 59 1 , the comoving spatial Laplacian of the Ray- 



(59) chaudhuri equation (38 1, the comoving spatial Laplacian 



of the Gauss-Codacci equation ( 45 ) and the spatial gra- 



dient of the momentum conservation equation, which to 
first-order reduce to, 

Ap + QA{p + p) + {p + p) AO = , (65) 

Ae + |eAe + |A(p + 3p)-A(£) 6 a fe ) =0 , (66) 

AJC + |6A6 - 2nAp = , (67) 

Ap + R 2 (p + p) D b a b = . (68) 

In order to reverse the order of the comoving spatial 
Laplacian and the time derivative of a scalar field S in 
the propagation equations ( 65 ) and ( 66 1 , it is useful to 



introduce the first-order scalar identity. 



AS = (AS)' 



±OAS 



R 2 SD b a b . 



(69) 



The dynamical equations ( 65 ) , ( 66 1 , and ( 67 1 are recast 
in terms of the comoving spatial Laplacian (611, (62 1 



and ( 63 1 using the the scalar identity (69 1, the gradient 



of the momentum conservation equation ( 68 1 , the Fried- 



mann equation ( 47 ) , the Raychaudhuri equation ( 38 1 and 



the energy conservation equation ( 59 1 . To first-order, in 



a comoving frame, the dynamical equations reduce to, 



+ np(l + w)$ = , 

1 3Kc 2 s 

2 + R 2 np(l + w) 



R?np(\ + w 



x+|e*-$ = o 



(70) 

-A$ = , 

(71) 
(72) 



In order to determine the dynamics of the scalar per- 
turbations, it is useful to express the dynamical equa- 



cording to, 



tions (70 1, (71 1 and (72 1 in terms of conformal time ac- 



+ H<Z> + Rnp(l + w)* = 
R 3Kc 2 s 
2 + Rkp(1 + w) , 



Rnp(l + w) 



2H 

x+— *-$ = 



(73) 

A$ = , 
(74) 

(75) 
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In order to determine the conformal time evolution of 
the density perturbation variable <f>, the dynamics has to 
be recast in terms of a second-order differential equation 
for By differentiating the ^-propagation equation ( 73 ) 



with respect to conformal time, using the ^-propagation 
equation (74 1 to substitute <& for if? and the zero-order re- 
lations (50 I and (51 1, a second-order differential equation 



for <f> is obtained according to, 



3H (1 + cf) $' 



K) + m') <f> - c 2 s A<f> = 



(76) 



which is the Bardeen equation, denoted [4.9] in [2] and 
identified by Woszczyna & Kulak [15 in the 1 + 3 covari- 
ant formalism. 

We now formally relate the energy density perturba- 
tion variable used in the 1+3 covariant approach to the 
Bardeen variables $a and <f># used in the 'background- 
based' approach. By taking the comoving spatial Lapla- 
cian of the divergence of the electric part of the Weyl 
tensor E a j, (56), the energy density perturbation <i> is 
found to be related to E a b 



as, 



$ = 3R 3 D a D b E„ 



(77) 



The expression of the electric part of the Weyl tensor 
E a f, in terms of the Bardeen variables $a and was 
first established by Bruni, Dunsby & Ellis [5] in equa- 
tions [113-114]. For a perfect fluid (i.e. in absence of 
anisotropic stress), the Bardeen variables have the same 
norm but opposite signs (i.e. &h = —&a)- Substitut- 
ing the relations [113-114] derived by Bruni, Dunsby & 
Ellis [5] into ( 77 1 , the energy density perturbation 4> is 
found to be the fourth-order derivative of the Bardeen 
variable <&a according to, 



where, 



□ = 3D a D D/ n D 



(78) 



(79) 



Note that by taking spatial derivatives of the electric part 
of the Weyl tensor, the vector and tensor perturbation 
terms contained in equations [113-114] vanish. To first- 
order, the conformal time derivative and the comoving 
spatial Laplacian commute with the spatial differential 
operator □, 



= n<s>' A , 

A$ = DA$i 



(80) 
(81) 



Thus, to first-order, the dynamics of the energy den- 
sity perturbation <I> is identical to the dynamics of the 
Bardeen variable <&a since the Bardeen equation (76) is 



entirely determined by conformal time derivatives and 
spatial Laplacians of <&. Therefore, in the 1 + 3 covariant 
formalism, $ is the analog of $a- 



B. Mukhanov-Sasaki equation 

The Bardeen equation describes the time evolution 
of the density perturbation variable Similarly, the 
Mukhanov-Sasaki equation describes the time evolution 
of the curvature perturbation variable, which in the 
spatially-curved case we will denote by £, and which 
reduces to \ defined in ( 63 1 in the spatially-flat case 
(K = 0), as we will now show. 

To make contact with the standard definition of the 
curvature perturbation in the 'background-based' ap- 
proach, it is useful to express \ m terms of $ only. By 
substituting ( 73 1 into ( 75 1 , using the Friedmann equa- 



tion (48 1, the expression for x to first-order becomes, 



X = 



3(1 + w) 



K 
W 2 



(82) 



which, for a spatially-flat background space-time (K 
0), reduces to, 



X 



3(1 +w) 



n 



(83) 



which is the standard definition of the curvature pertur- 
bation (see Mukhanov [57], Durrer [2H])- In the spatially- 
flat case, the conformal time derivative of the curvature 
perturbation \ is given by, 



X = 



2cg 
3W(1 + 



w ) 



A$ 



(84) 



which shows that the curvature perturbation x is a con- 
served quantity on large scales. To show this explicitly a 
Fourier transform has to be performed (see Durrer [28). 
However, we note that on large scales the comoving spa- 
tial Laplacian of the analog of the Bardeen variable is 
negligible compared to the comoving scale TL, thus satis- 
fying A<I> <C Tt- Nonetheless, the curvature perturbation 
is not conserved on super-Hubble scales for a spatially- 
curved background space-time [5]. 

In the spatially-curved case, Bruni, Dunsby & Ellis [5j 
mention that there is a generalised curvature perturba- 
tion C, which is conserved on large scales and defined up 
to a constant amplitude. Thus, we define the generalised 
curvature variable £ as 



r- R f A! r 4K A P 
^ ~ 2 V R 2 (l + w) P 



(85) 



where £ = \C in order for £ to reduce to x m the 
spatially-flat case (K = 0). The generalised curvature 
perturbation variable £ in the spatially-curved case is re- 
lated to the curvature perturbation \ i n the spatially-flat 
case by, 



C = x- 



2K 



3H 2 (l 



K 



(86) 
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We now briefly show that £ is conserved on large scales. 
By substituting ( 82 1 into ( 86 1 , the generalised curvature 



perturbation £ is recast in terms of the analog of the 
Bardeen variable <i> only according to, 



c 



3(1 



K 
V 2 



K 



n 



For a spatially-curved case, by differentiating (87) and 



substituting the Bardeen equation ( 76 ) , the conformal 



time derivative of the generalised curvature perturbation 
£' is found to be, 



2c? 



m{i+w) 



K 

n 2 



which means that the generalised curvature perturba- 
tion £ is a conserved quantity on super-Hubble scales. 
To show this explicitly a harmonic decomposition has to 
be performed. However, similarly to the spatially-flat 
case, we note that on large scale the comoving spacial 
Laplacian of the analog of the Bardeen variable is negli- 
gible compared to the comoving scale TL, thus satisfying 
A$ < TL. 

In order to determine the conformal time evolution of 
the curvature perturbation variable £, the dynamics has 
to be recast in terms of a second-order differential equa- 
tion for Q. By differentiating the £ _ P r opagation equa- 
tion ( 88 1 with respect to time and using the comov- 



ing Laplacian of the generalised curvature perturbation 
C (87 1 to substitute ( for <f>, a second order differential 



equation for £ is obtained according to, 



->!?(> 



3d 



TL 



C - c 2 s A( = 



It is convenient to introduce a new variable z, 



R 2 (K 



c s TL 



1/2 



(89) 



(90) 



which allows the ^-propagation equation ( 89 1 to be recast 
in a more convenient form as, 



C" + 2-c'- c 2ac = o 



(91) 



Note that using the dynamical equation ( 50 1 , z is rewrit- 
ten as, 



R?_ 

c s n 



(H 2 -H' + K 



1/2 



(92) 



which corresponds to the variable defined by Mukhanov 
et al. [5] in equation [10.43b]. 

It is also useful to define another variable, 



v = zC, 



(93) 



which is the analog of the variable v defined by Mukhanov 
et al. [3] in equation [10.61]. Using v, the dynamics of the 



second-order dynamical equation (91 ) is explicitly recast 



in terms of a wave equation given by, 







(94) 



This wave equation is the Mukhanov-Sasaki equation in a 
spatially-curved FLRW background space-time denoted 
[11.7] in 0. 



C. Mukhanov-Sasaki variable 

The scalar variable v = z( is the analog of the 
Mukhanov-Sasaki variable associated to the scalar per- 
turbations [21 3] (note that in a spatially-flat case, the 
analog of the Mukhanov-Sasaki variable reduces to v — 
zx). In the general case, the homogeneity and isotropy 
of the spatial hypersurfaces enable us to perform a har- 
monic decomposition of the scalar perturbation variable 
v, where v is decomposed into components that trans- 
form irreducibly under translations and rotations, and 
evolve independently, as explained in [28 . The harmonic 
analysis of the scalar perturbation v consists of a de- 
composition into eigenfunctions of the comoving spatial 
Laplacian of v according to, 



-k 2 



Vk 



(95) 



where k is the eigenvalue of the associated harmonic 
mode and the fc-index denotes the eigenvector of the 
mode. The comoving wavenumber v of the scalar mode 
is defined as, 



v 2 = k 2 + K , 



(96) 



where K = { — 1,0,1} is normalised. The comov- 
ing wavenumber v takes continuous values when K — 
{ — 1,0} and discrete ones for K = 1. In particular, 
the regular normalisable eigenmodes have v > for 
flat and hyperbolic spatial hypersurfaces, and an integer 
satisfying v > 1 for spheric hypersurfaces as explained 
in [23l|29]. 

To first order, the dynamics of the scalar perturba- 
tions can be rewritten as a serie of decoupled harmonic 
oscillators. Using the harmonic decomposition (95 1 in 



terms of k, the Mukhanov-Sasaki wave equation |94|) in 
the A;-mode is given by, 



:k 



z 
z 



v k = . 



(97) 



Note that ( 97 ) corresponds to a simple harmonic oscilla- 



tor in the fc-mode, 



v k " + u\v k = , 



(98) 



where the conformal time dependent frequency Wfc(f) is 
given by, 



LO k = c s k 



,2 ;,2 



1/2 



(99) 
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D. Comparison with scalar perturbations of a 
scalar field in the 'background-based' approach 

The dynamics of the early universe is believed to un- 
dergo an inflation phase described by a scalar field [3UJ 
I5T] . In a 'background-based' approach, the scalar field <f> 
in the 'real' space-time is decomposed into a background 
component <po and a gauge-dependent perturbation 5<p 
according to, 



0(i,x) = 0o(t) + <ty(*,x) • 



(100) 



According to Weinberg [32j . the energy-momentum ten- 
sor of an unperturbed scalar field takes the perfect fluid 
form with an energy density and pressure respectively 
given by, 



Po 
Po 



¥1 



V(<fo) 



(101) 



where V(4>o) is an arbitrary real potential. The initial 
conditions for inflation needed to perform the quantisa- 
tion are set as the model emerges from the Big Bang, 
where the scalar field dynamics is dominated by the ki- 
netic term and satisfies <f>Q 3> V(c/)q). Therefore, at very 
early times, the scalar field behaves like stiff matter and 
the equation of state parameter w is related to the speed 
of sound by w — c 2 = 1 [28] . The dynamics of scalar per- 
turbations of a scalar field in a 1 + 3 covariant approach 
was first investigated by Bruni, Ellis & Dunsby [33 . Lan- 
glois& Vernizzi [53] later generalised the 1 + 3 covariant 
approach to multi-scalar fields. 

To first order, the dynamical equation for the Bardeen 
variable is given by, 



(102) 



which corresponds to equation [10.1.12] presented by 
Weinberg [32 . By comparing this dynamical relation 
with the corresponding result ( 73 1 found in the 1 + 3 co- 



variant approach using the zero-order expression for the 
energy density and pressure of a scalar field ( 101 1 and by 
using the correspondence relation ( 78 1 , the 1 + 3 covari- 



ant perturbation variable 'J is found to be related to the 
scalar field perturbation 5<fi by, 



-□ 



20 o 



(103) 



The curvature perturbation variables in the flat case x 
and the curved case £ can be recast in terms of the 
Bardeen variable <j>^ and the scalar field perturbation 
variable 5(f> according to, 



x = n [$a + h-^ 



C = a 



u 

K f K 

m 2 { + n 2 



(104) 



The quantisation variable z can also be expressed in 
terms of the scalar field perturbation according to, 



(!) 



(105) 



Finally, for a scalar field behaving like a stiff fluid (i.e. 
c s = 1), the Mukhanov-Sasaki wave equation in the k- 
mode (97 1 reduces to, 



Vk 



(106) 



Thereby, with the identities presented above, we formally 
related the first-order scalar perturbation variables in the 
1+3 covariant formalism to the corresponding variables 
in the 'background-based' approach. 



V. VECTOR PERTURBATIONS 

Vector perturbations are described by spatially pro- 
jected and divergence- free vectors |29j . The only dynam- 
ical variable which satisfies these constraints to first-order 
is the vorticity pseudo- vector uo a . The vanishing diver- 
gence of uj a to first-order can be deduced from (41 1. 

The dynamics of the vorticity covector uj a is deter- 
mined by the vorticity propagation equation ( 39 1 . To 



obtain an explicit evolution equation in terms of the vor- 
ticity covector only, the term involving the acceleration 
has to be recast in terms of the vorticity. Using the mo- 
mentum conservation equation ( |60[ ) and the kinetic non- 
commutation identity (13), we find to first-order, 



curla„ = 2c o 0t/ 



(107) 



Thus, to first-order, the vorticity propagation equa- 
tion (39 1 can be recast as, 



uj (a) + |6 (1 - fcf) u a = , (108) 

as shown by Hawking [10J. It is convenient to express 
this dynamical equation in a comoving frame according 
to, 



J (a) 



m (i 



3„2 



<S Wa = 



Using the vorticity contraction identity, 



w 2 ^2Lo[ a) uj a 



(109) 



(110) 



a first-order propagation equation for the vorticity scalar 
w is obtained. 



= -2H(l- %ci)w , 



(HI) 

and implies that the vorticity scalar scales as [TOl , 

w «r 2+3c -. (ii2) 
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The scaling relation (1121 implies that during the ex- 



pansion phase (R' > 0), the scalar amplitude of the vor- 
ticity u decays if < |. The inflationary scenario is 
the simplest known generating mechanism for the ini- 
tial density fluctuations. For an inflaton field in slow-roll 
(w = —1), the vorticity scalar scales as ui <x R~ 5 . Hence, 
the value of the vorticity scalar at the end of slow-roll 
inflation wt is related to the value of the vorticity scalar 
at the onset of the slow-roll inflation u)i by, 



exp (-5JV) 



(113) 



where N is the number of e-fold during the slow-roll 
phase. Hence, if the vector perturbations were initially 
significant, they have decayed by a factor exp(5iV) dur- 
ing slow-roll inflation (113) and can safely be neglected 



in a subsequent quantitative perturbation analysis. 



VI. TENSOR PERTURBATIONS 

Tensor perturbations are described by spatially pro- 
jected, symmetric, trace-free and transverse second rank 
tensors, as explained in To find a suitable tensor to 
describe such perturbations, it is useful to split the mag- 
netic part of the Weyl tensor H ab into a transverse part 



denoted H^P and a non-transverse part denoted H K a 
according to, 



(V) 



H 



(T) 



H 



(V) 



(114) 



where the ^ and indices refer respectively to tenso- 
rial and vectorial degrees of freedom. By construction, 

(T) 

Kb 

is divergence-free and satisfies the requirements for 
a tensor perturbation. For convenience, in this section, 
we will not use explicitly the ( - T - ) index to refer to the 
transverse part of the magnetic part of the Weyl tensor, 
since we are only considering tensorial degrees of free- 
dom. Note that in the irrotational case, the magnetic 
part of the Weyl tensor is divergence-free to first-order 
— this can be deduced from the constraint (57 1. Thus, 



in that case, the analysis is restricted to the tensorial de- 

(T) 

grees of freedom only and H ab = H y ab ' . It is also worth 
mentioning that the electric part of the Weyl tensor E ab 
cannot qualify as a tensor perturbation, since it is not 
divergence-free in presence of matter — this can be in- 
ferred from the constraint ( 56 1 . 



A. Grishchuk equation 

In order to obtain a second-order differential equation 
in terms of H ab , it is useful to introduce linearised iden- 
tities found in [13] • Using the Ricci identities (36l-(37l, 



the expression for the Riemann tensor to zero-order ( 34 ) 
and the definition of a curl (10 1, a symmetric and spa- 
tially projected tensor T/ a b) to first order has to satisfy 



the geometric linearised identity, 

(3) (curlT ab )' = curlT (ah) - ±6curlT afc 
and if T ab is transverse to first-order, 

D b T ab = , 
then it also satisfies the identity, 



curl ( curl T ab ) = -D 2 T ab 



3K r 

i?2 



(115) 
(116) 

(117) 



By differentiating the magnetic propagation equa- 
tion ( 55 1 , substituting the electric propagation equa- 



tion (54 and using the linearised identities (1151 



and (117), the dynamical equation for H ab to first or- 



der reduces to, 

7 / 3K\ 

H {ab) + -QH {ab) +2 ( (1 - w)k P - - H ab -D 2 H ab = , 

(118) 

which has been established in [TS] . Using the zero-order 
dynamical relations ( 48 ) and ( 49 ) , it is convenient to re- 



cast the dynamics of H ab (1181 in terms of conformal 
time f according to, 



u" 



6HH', 



(ab) 



(A — 2K — m 2 - AH 1 ) H ab = . 

(119) 



In order to eliminate the second term on the LHS of ( |119 1 
and find a suitable wave equation for the gravitational 
waves, it is useful to introduce the rescaled magnetic part 
of the Weyl tensor H ab defined as 



R H ab 



(120) 



The propagation equation of the tensor perturba- 
tions (1191 can be now elegantly reformulated to first- 



order in terms of H ab and yields 

R" 



11 (ab) 



2K 



R 



, 



(121) 



which is our version of the Grishchuk equation |20j de- 
scribing the dynamics of primordial gravitational waves 
in a spatially-curved case. 



B. Tensor decomposition of H ab 

In order to determine the scalar canonical variables as- 
sociated to the tensor perturbations, which are the scalar 
amplitudes of the tensor perturbations and will be re- 
ferred to as canonical variables |20j in the current work, a 
tensor decomposition of H ab has to be performed. From 
the properties of the magnetic part of the Weyl tensor 
H abl we deduce that the key tensor H ab is symmetric, 
trace-free and transverse to first-order. The transversal- 
ity of H ab can be expressed as, 



D b H n 



k b H„ 







(122) 
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where k b is the spatial wavevector of the transverse grav- 
itational waves satisfying u b kb — 0. For convenience, we 
now define two vectors e^ a and e^ a that provide an 
orthonormal basis for the two-dimensional spatial hyper- 
surface orthogonal to the propagation direction k a of the 
gravitational waves, and thus satisfy the following con- 
straints, 



fc«e« = rei 2 > = , 
P (i)<vW - p( 2 ) a p( 2 ) - 1 

e« a e ( 2 > = . 



(123) 



Taking the comoving Laplacian of the orthonormality 
conditions ( | 128 1 and using the identities (1261, we de- 
duce the following constraints, 



e+Ae*<* 



e„,Ae 

ab 



lib 



Xab 

Ae +ab 



, 
. 



(130) 



The tensor perturbation variable H a b is spatially pro- 
jected, transverse, traceless, and therefore can be decom- 
posed into two polarisation modes {+, x} according to, 



(131) 



where 



Note that these two vectors are not uniquely defined, 
which does not hinder our perturbation analysis since 
any orthonormal vector basis of the two dimensional hy- 
persurface can be rotated to recover our vector basis 

| e (l)a je (2) a |_ 

For an irrotational fluid, the Fermi- Walker transport 
of the basis vectors el and e a 2%> vanishes 1351 . 



u a a b e ( ^ = 



u b V b eW-u a a b eW = n 



(124) 



Projecting the Fermi-Walker transported basis vec- 
tors (1241 on the spatial hypersurface yields, 



h x = H a 



j^jab c + 
b„x 



b • 



J ab ' 



are the scalar amplitudes of the tensor perturbations. 
The corresponding decomposition in the 'background- 
based' approach is mentioned by Durrer in |28j . The lin- 
ear independence of the polarisation tensors in the first- 
order perturbation analysis allows us to study separately 
the dynamics of the two decoupled polarisation modes. 
To keep the notation compact, we introduce the polari- 
sation mode superscript A = {+, x}. 



Canonical variables 



(iy _ (2)/ _ 

e (a) ~ e (a) - U 



(125) 



Taking the comoving spatial Laplacian of the con- 
straints (1231, we deduce the following identities, 



e d)a Ae (l) = g (2)a Ae (2) = q ^ 
e (l)a Ae (2) = _ e (2)a Ae (l) _ 



(126) 



In order to decompose the tensor perturbation tensor 
H a b into two polarisation modes, we define two covariant, 
trace-free and linearly independent polarisation tensors, 



26 2 
1 



1 LWpW _ P (2)J2) 



'a °6 



, - _ , e (l) e ( 2 ) 4- e ( 2 ) e (1) 
C ab ~ 2 I e a C b T e a e & 



(127) 



which satisfy the orthonormality conditions, 



P . P , ab _ X xab _ i 



-ab z 
Z ab e 



e\e xab = 



(128) 



Contracting the Grishchuk equation ( 121 1 with the po- 



larisation basis tensors e\ h and using the identities (jl 29 1 



and (1301, we obtain the Grishchuk equation for the 



scalar amplitude of the tensor perturbations associated 
to the polarisation mode A, 



Xll 



A-2K 



IT 

If 







(132) 



where h x are the canonical variables, which describe the 
scalar amplitude of the tensor perturbations in the po- 
larisation modes A. There are two canonical variables 
{h + ,h x } associated to the polarisation modes of the ten- 
sor perturbations. 

In a similar manner than for the scalar perturbations, 
we perform a harmonic decomposition of the tensors per- 
turbations. The harmonics analysis of the tensor pertur- 
bations consists in a decomposition into eigenfunctions of 
the comoving spatial Laplacian of the scalar amplitude 
of the tensor perturbation variable h x according to, 



2. A 



-k Z h 



(133) 



and thus form an orthonormal tensor basis for the polar- 
isation modes From the identities (125), we deduce 



that the spatially projected time derivatives of the polar- 
isation tensors vanish, 



-{ab) 



ef' = 

(ab) 



(129) 



where k is the eigenvalue of the associated harmonic 
mode and the suffix k denotes the eigenvector of the 
mode. 

Finally, to first-order, the dynamics of the tensor per- 
turbations for each polarisation mode A can be rewrit- 
ten as a serie of decoupled harmonic oscillators. Using 
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the harmonic decomposition of the scalar amplitudes in 



terms of k ( 133 1 , the evolution equation for the scalar am- 



plitude of the tensor perturbations ( 132 1 in the fc-mode 
satisfies, 



k 2 + 2K 



R 



(134) 



for each polarisation A = {+, x}. The tensor pertur- 
bation amplitudes variables h x are identified as the two 
scalar canonical variables associated to the tensor per- 
turbations [20] . Note that ( 134 1 corresponds to a simple 



harmonic oscillator with a polarisation A in the fc-mode, 



■ LV 



o , 



(135) 



where the conformal time dependent frequency <x>£(t) is 
given by, 



2K 



1/2 



(136) 



the vector perturbations, which does not lead to a sec- 
ond order wave-equation, unlike the scalar and tensor 
perturbation, but to a first order equation involving the 
vorticity. In an expanding universe, the vorticity decays 
provided the sound speed squared satisfies c 2 < | and 
the vector perturbations can therefore be neglected in 
the analysis. 

In cosmological perturbation theory, there are in gen- 
eral six degrees of freedom: two scalars, two vectors and 
two tensors as mentioned by Bertschinger |36j . How- 
ever, in the perfect fluid case, the two Bardeen potential 
are related by <&a = — $h, and thus there is only one 
remaining scalar degree of freedom. To sum up, to first- 
order, there is one degree of freedom associated with the 
scalar analog of the Mukhanov-Sasaki variable v, two de- 
grees of freedom associated with the divergence-free vor- 
ticity covector u> a and two degrees of freedom related to 
the two canonical variables h + and h x representing the 
scalar amplitudes of the tensor perturbations. Hence, to 
first-order, the dynamics of the adiabatic perturbations 
of a perfect fluid are described by five parameters, as 
expected. 



VII. CONCLUSION 



We performed a perturbation analysis of an adiabatic 
perfect fluid to first order using the 1 + 3 covariant and 
gauge-invariant formalism and identified the analog of 
the Mukhanov-Sasaki variable and the canonical vari- 
ables needed to quantise the scalar and tensor perturba- 
tions respectively about a spatially-curved FLRW back- 
ground space-time. We also determined the dynamics of 
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